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Millimeter observations of obscured Spitzer 24 micron sources 
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ABSTRACT 

We present MAMBO 1.2mm observations of 40 extragalactic sources from the 
Spitzer First Look Survey that are bright in the mid-infrared (S 2 4 Mm > lmJy) 
but optically obscured (log\Q(vF v (2Aixrri) / vF v {Q.7 ixm)) > 1). We use these ob- 
servations to search for cold dust emission, probing the similarity of their spec- 
tral energy distributions to star forming infrared galaxies or obscured AGN. 
The sample as a whole is well detected at mean Si. 2mm = 0.74 ± 0.09mJy and 
Si.2mm/S24 A tm = 0.15 ± 0.03. Seven (three) of the sources are individually de- 
tected at > 3(7 (> 5cr) levels. Mean millimeter fluxes are higher for sources with 
the reddest mid-infrared/optical colors. Optically faint but with relatively low 
mm to mid-infrared ratio, the typical SEDs are inconsistent with redshifted SED 
shapes of local star-forming infrared galaxies. They also differ from SEDs of typ- 
ical submillimeter selected galaxies, with the 24/im sources that are individually 
detected by MAMBO possibly representing intermediate objects. Compared to 
star-forming galaxies, a stronger but optically obscured mid-infrared component 
without associated strong far-infrared emission has to be included. This compo- 
nent may be due to luminous optically obscured AGN, which would represent a 
significant part of the high redshift AGN population. 
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1. Introduction 

Significant progress has been made in the last few years in detecting infrared galaxies 
at high redshift and characterizing their nature, evolution, and contribution to the cosmic 
background. ISOCAM has detected at 15 fim a population of z < 1 luminous star-forming 
galaxies (e.g. Elbaz et al. 2002), still fairly bright in the optical (e.g. Flores et al. 1999). 
SCUBA and MAMBO surveys of small fields detect a distinct z ~ 2.5 population of hy- 
perluminous star-forming and active galaxies (for a review see Blain et al. 2002) that are 
optically extremely faint (e.g. Ivison et al. 2002; Dannerbauer et al. 2004). Spitzer Space 
Telescope 24 /im surveys sample the redshifted, rest-frame 6 — 12 /xm strong dust emission 
from PAH and very small grains in galaxies at z ~ 1 — 3. Therefore, deep mid-IR observations 
using MIPS on Spitzer offer a new and efficient probe of dusty galaxies at high redshift (e.g. 
Perez Gonzalez et al. 2005). Spitzer 24/zm surveys are expected to detect all of the ISOCAM 
and most of the SCUBA/MAMBO sources, bridge the gap between these two populations 
in redshift and luminosity, and search for entirely new categories of sources. 

The first non-proprietary Spitzer survey is the First Look Survey (FLS), covering 3.7 
square degrees in seven bands. More than 18000 sources are detected to the 3a sensitivity of 
110 /xJy at 24 /im. Combining deep optical, Spitzer IRAC 8 /im, and MIPS 24 jum data, Yan et 
al. (2004) have obtained the first characterization of a large sample of 24 /im-selected sources. 
These sources have colors ranging from those of relatively low-redshift, fairly inactive galaxies 
to extremely red 24/xm/i? and 24/im/8/xm colors that can only be explained by luminous 
obscured starbursts or AGN at high redshift. The available color information is not unique 
for identifying the nature of these objects, however, and the applicability of the local SED 
templates has to be critically examined. This is particularly true for the rare and extreme 
24 /im-selected objects that the FLS sample includes due to its large size. In this letter we 
report on MAMBO 1.2mm observations of such rare 24/im-bright but optically faint sources, 
providing a constraint on the cold dust contribution to their spectral energy distributions 
that is essential for characterizing their nature and their relation to other populations of 
high-redshift infrared galaxies. 



2. Sample and Results 

Yan et al. (2004) used Spitzer mid-IR and ground-based optical photometry for a first 
characterization of 24/mi sources from the FLS. We have selected for 1.2mm photometry a 
subset of bright (S 2 4 Mm , >lmJy) sources that are optically faint, as evidenced by a 24/zm 
to optical 'color' #(24, 0.7) = log w {vF v {2^m)/vF v {Q.l //m)) of at least 1 (R>22.5 for 
S24/im =lmJy). We have chosen 39 of the 145 sources meeting these criteria, and added 



- 3- 



one with a slightly lower R(24,0.7). We refer to these sources as 'obscured 24/im sources'. 
The selection covers the full color spread of the Yan et al. (2004) objects (Fig. 1), preferring 
objects brightest at 24/im, and aiming for a large overlap with FLS objects studied by Spitzer 
mid-IR spectroscopy (Yan et al. 2005, and in preparation). Objects were selected using the 
Yan et al. (2004) flux catalogs. Since then, slightly revised 24/im and 8/im fluxes (M. Lacy 
et al. in preparation) have become available which we adopt in the following discussion. 

Observations were spread over the pool observing sessions at the IRAM 30m telescope in 
fall/winter 2004/2005, using the 117 element version of the Max Planck Millimeter Bolometer 
(MAMBO) array (Kreysa et al. 1998) operating at a wavelength of 1.2mm. On-off observa- 
tions were typically obtained in blocks of 6 scans of 20 subscans each, and repeated in later 
observing nights until reaching noise near 0.66mJy or a 5a detection. The median 1.2mm 
noise level reached for the sample is 0.59mJy (range 0.35 to 0.82mJy). The data were re- 
duced with standard procedures in the MOPSIC package developed by R. Zylka, using the 
default calibration files for these two pool periods. Table 1 lists the resulting 1.2mm fluxes 
and their statistical uncertainties. 

Three sources are detected individually at more than 5a with 1.2mm fluxes between 
2.3 and 5.8mJy, and four more at 3a to 5a. The total 40 object sample is well detected 
statistically at a mean flux of 0.74±0.09mJy (Table 2). To investigate a dependence on the 
mid-IR/optical colors, we have grouped the sample in four color 'quadrants' separated at 
R(24,0.7)=1.5 and and R(24,8)=0.5. Due to the sample selection, these quadrants contain 
similar numbers of objects by design. Table 2 shows clear trends in the mean 1.2mm fluxes 
for these subsamples. Objects redder in R(24,0.7) have significantly higher mean 1.2mm flux, 
and a trend may also be present with R(24,8). Values quoted in Table 2 refer to means of 
the individual source data weighted by the inverse of the noise squared. Very similar values 
are obtained for equal weight means. The trend in the mean fluxes agrees with the > 3a 
detection of 7/24 objects at R(24,0.7)>1.5 vs. 0/16 at bluer color. Using the formalism of 
Stevens et al. (2005), this implies that the two groups are different at 2.9<r significance. The 
distribution of 1.2mm vs. 24/im fluxes, shown in Fig. 2, is inconsistent with a single flux 
ratio for the entire sample, or for the subsample with R(24,0.7)>1.5. This is expected given 
the sample selection, and emphasizes that the mean values quoted in Table 2 are formed 
over (sub) samples that can have noticeable spreads of intrinsic properties. 

3. Millimeter properties of obscured Spitzer sources 

Our color selection implies that these objects must be heavily obscured and/or at con- 
siderable redshift. This is illustrated in Fig. 3, showing that even spectral energy distribu- 
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tions similar to 'standard' dusty starbursts like M82 fail to reproduce for plausible redshifts 
R(24,0.7)>1.2, i.e. the vast majority of our sample. R(24,0.7) will be yet lower for SEDs of 
UV-brighter starbursts like NGC 7714. SEDs similar to the most optically obscured local 
ultraluminous infrared galaxies like Arp220 reach R(24,0.7)>1.5 for redshifts around 2, but 
even stronger emission at observed 24/xm (rest-frame 8/xm for z«2) relative to the shorter 
wavelengths is needed to reproduce the most extreme Spitzer objects. The same is true for 
the SED of the well studied less luminous but extremely obscured galaxy NGC 4418. For 
this object, both the 8/zm and the 0.7/xm flux are host dominated at z>1.5. Varying the 
nucleus to host ratio moves the source diagonally in Fig. 3. 

The millimeter observations provide a crucial test of the analogy to local SEDs. Fig. 4 
shows the expected ratio of 1.2mm and 24/iin flux density as a function of redshift, again for 
M82, Arp220, and NGC 4418 SEDs. The seven objects with individual MAMBO detections 
are reasonably close to the properties of these local templates. At S\.2mml Si^m ~ 1 . . .5, 
they are in a range still consistent with local infrared galaxy SEDs redshifted to z~2, assum- 
ing a relatively high optical obscuration or modest additional mid-iR emission. In contrast, 
comparison of the templates to the mean values for the obscured Spitzer sources (Table 2) 
shows that in fact neither SED provides a good representation of the mean properties of our 
24/im sources. Arp 220 and NGC 4418-like SEDs, which were in a reasonable range of the 
optical/mi d-IR colors, predict by far too strong 1.2mm flux for plausible redshifts. M82-like 
SEDs are closer to the observed mm to MIR ratio, but predict too bright optical counter- 
parts as noted above. We conclude from the analysis of mean properties that the sample of 
obscured 24/zm contains many sources that are not well fit by these local infrared luminous 
galaxy SED templates. Relative to the templates, they require additional rest frame mid-IR 
emission, due to a component that is both warm (to avoid overproducing mm emission) and 
optically obscured. 

Another way of phrasing the need for 'additional rest frame mid-IR emission' is to state 
that the overall infrared SED must be warm to hot, without a dominant rest frame far-IR 
peak. In a strongly simplifying single temperature blackbody fit, the observed 1.2mm to 
24/im flux density ratio corresponds to a temperature of ~75K, i.e. intrinsic T>200K if at 
z~2, and luminosities of the order of 10 13 L Q . The warm temperature is suggestive of a strong 
AGN contribution. The SED of the most luminous radio-quiet local quasar, PDS 456 (Yun 
et al. 2004) is able to fit the 1.2mm/24/im ratio for plausible redshifts (Fig. 4), but this Type 
1 object is too bright in the UV/optical to explain the mid-IR/optical colors for our sample. 
That more obscured AGN may be able to meet that constraint can be illustrated using the 
nuclear SED of the classical Seyfert 2 NGC 1068. Its nuclear FIR/submm emission can 
currently be separated from the host only at 450yum with 8450^™ ~ l-5Jy (Papadopoulos & 
Seaquist 1999), but this can be combined with the mid-IR AGN continuum (Lutz et al. 2000) 
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and the rest frame UV taken from NED, to conclude that such an AGN SED would satisfy 
the criteria for both faintness at 1.2mm (Fig. 4) and optical obscuration (Fig. 1). AGN in 
obscured 24/xm sources must suffer intermediate obscuration, sufficient to suppress the rest 
frame UV/optical but not moving too much energy into far-IR re-emission. Separation from 
host emission and spectral characterization of the rest frame far-IR emission is currently 
incomplete for local obscured AGN, but objects like IRAS F00183-00183 (Tran et al. 2001; 
Spoon et al. 2004) do have more favourable mid- to far-IR ratios than NGC 4418 and may 
fit the SEDs of some obscured 24/im sources. 

Low 1.2mm/24/im flux ratios can also be obtained for low redshift (z<l) infrared star- 
forming galaxy SEDs. As already noted, this is ruled out for the adopted templates, by the 
failure to meet R(24,0.7)>1 at such redshifts (Fig. 3). Contamination by similar sources 
may occur at the low R(24,0.7) end of the MAMBO sample if the scatter in rest frame SEDs 
extends towards smaller optical/mi d-IR flux ratios. This scenario becomes progressively 
unlikely for the R(24,0.7)>1.5 sources with weak 1.2mm emission. 

Hunt et al. (2005) use the likely compactness of star formation events in high redshift 
infrared galaxies to argue for SEDs of Blue Compact Dwarf Galaxies as local analogues. The 
most extreme case, SBS 0335-052, shows an absorbed warm mid-IR continuum (Thuan et al. 
1999; Houck et al. 2004). While this possible analogy has to be mentioned, it is unclear to 
which extent such low metallicity dwarfs can relate to the luminous obscured 24/im sources. 
SBS 0335-052, in addition to its mid-IR emission, has strong unobscured optical/UV emission 
which would place it far off the color selection of obscured 24/im sources. Rest frame submm 
information is not available to place SBS 0335-052 on Fig. 4. 

4. Relation to other populations 

(Sub)millimeter galaxies (SMGs) are the best studied pre-Spitzer population of z> 1.5 
infrared galaxies, inferred to be objects at z~2.5 with roughly ULIRG like SEDs but much 
higher luminosity. Optical to mid-IR colors of bright and obscured 24/im sources (Figs. 1,3) 
and their low surface density compared to that of SMGs could have been consistent with an 
identification of the yet more extreme luminosity end of the SMG population as obscured 
24/im sources. The small ratio of 1.2mm and 24/im emission clearly argues against this being 
true for the majority of the Spitzer sources. Egami et al. (2004), Ivison et al. (2004), and 
Frayer et al. (2004) have reported 24/mi photometry of SMGs from SCUBA and MAMBO 
blank field and radio preselected surveys. They find much higher Si,2 mm / 5*24^™, with median 
~ 12, clearly inconsistent with SMGs' and obscured Spitzer sources' having similar SED 
shapes and redshifts. Again, the Spitzer objects with individual MAMBO detections are 
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closer to the SMG properties. 

Blain et al. (2004) discuss selection effects against warm SEDs in the SMG population 
and the potential of Spitzer in avoiding them. In a related effort, Chapman et al. (2004) 
have obtained spectroscopic redshifts of optically faint /zJy radio sources without submm 
detections, inferring under the assumption of the radio - far-IR correlation that these must 
represent a population of ultraluminous star-forming z~2.2 galaxies characterized by intrin- 
sically warm or hot dust continua. Little is known at present about the rest-frame mid-IR 
properties of these sources, making an assessment of the overlap with obscured 24/zm sources 
difficult. Frayer et al. (2004) report faint 24/zm fluxes or nondetections of several optically 
faint radio sources that have SCUBA 850/zm flux upper limits of order 8mJy. If these results 
are representative for the full population of optically faint radio sources, then their possible 
warm FIR SEDs cannot extend down into the mid-IR, and their nature must be different 
from the obscured 24/im sources. If, for example, their radio emission is AGN boosted, their 
dust emission may be comparatively weak at all wavelengths. 

Mid-IR spectroscopy will play a key role in identifying the nature of the rare obscured 
24/im sources and their strong mid-IR emission. Houck et al. (2005) and Yan et al. (2005) 
have reported the first Spitzer IRS spectroscopy of such sources, finding a predominance of 
continua and absorbed continua and indications for starburst PAH emission in a subset of 
sources. The obscured continuum spectra, which are seen only in a small minority of local 
infrared galaxies (Spoon et al. 2002) but a third to half of the obscured 24/im sources in the 
subsamples already published by Houck et al. (2005) and Yan et al. (2005), indicate heavily 
embedded energy sources that might be AGN. Redshift estimates, where they can be derived, 
suggest z~1.5-3. As for the photometric evidence, comparison to initial Spitzer spectroscopy 
of SMGs which shows PAHs (Lutz et al. 2005) suggests that the two populations have only 
modest overlap, and that there is a large contribution of warm and possibly AGN like sources 
to the obscured 24/im population. 

If a good part of the bright obscured 24/im sources contain luminous high z AGN they 
will represent a significant component of the high redshift AGN population. Our MAMBO 
sample is a subset of about 150 similar objects (Fig. 1) detected over the 3.7 square degree 
FLS field. If many of these objects are Lb i ~ 1O 13 L AGN at z~1.5-3 their comoving number 
density will be of the order 10~ 6 Mpc~ 3 , comparable to that of similar bolometric luminosity 
(Mb < —25) optically selected quasars (Croom et al. 2004). The ongoing characterization 
of their redshift distribution and SEDs will shed more light on this issue. 

Based on observations carried out with the IRAM 30m telescope and with the Spitzer 
Space Telescope. We thank the IRAM staff for support and advice, all pool observers for 
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their help in obtaining the observations, and the referee for helpful comments. 
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Table 1. MAMBO 1.2mm fluxes of obscured 24/zm sources 
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Table 1 — Continued 
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17:20:05.96 


59:17:45.0 


0.39±0.72 


1.59 


22404 


17:21:51.77 


58:53:27.7 


0.66±0.59 


1.57 


22563 


17:22:49.00 


58:50:27.5 


5.79±0.82 


1.14 


22582 


17:21:24.58 


59:20:29.5 


2.22±0.53 


1.11 
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Table 2. Mean 1.2mm properties for different ranges of mid-IR/optical colors 



R(24,0.7) R(24,8) N a x . 2 Si. 2mm S 1 . 2mm /S 24 , tm 

mJy mJy 

all all 40 0.59 0.74±0.09 0.15±0.03 

<1.5 <0.5 8 0.65 0.15±0.22 0.01±0.05 

<1.5 >0.5 8 0.59 0.46±0.20 0.26±0.11 

>1.5 <0.5 15 0.57 0.91±0.14 0.18±0.05 

>1.5 >0.5 9 0.53 1.08±0.18 0.52±0.10 



Note. — See Fig. 1 for the distribution of sources within the 
R(24,0.7) and R(24,8) colors. <7i.2 indicates the median 1.2mm 
flux error for the sources in each bin. 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 

R(24,0.7) = log 10 (i/f„(24 M m)/i/f„(0.7 M m)) 

Fig. 1. — Optical/mi d-IR color-color plot of bright (>lmJy) 24/xm sources from the Yan et 
al. (2004) Spitzer FLS sample, indicated by small crosses. Sources with R(24,0.7)<1 are not 
shown. Asterisks highlight the targets observed by MAMBO, big asterisks those individually 
detected by MAMBO at the > 3a level. 8/xm and optical limits are 3cr. 



a: 



-0.5 



13 




2 4 6 8 10 

24 J iim flux density [mjy] 



Fig. 2. — MAMBO 1.2mm fluxes as a function of 24/im fluxes. 24/zm errors are smaller than 
the symbols. Sources with R(24,0.7)>1.5 are shown as diamonds, others as asterisks. 
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Fig. 3. — Mid-IR/optical colors for redshifted objects with the SEDs of the star forming 
infrared galaxies M82 and Arp220, and of the heavily obscured galaxy NGC 4418. Locations 
of the FLS sources observed by MAMBO are repeated from Fig. 1. Redshifts from to 3 are 
marked in intervals of 0.5. The adopted SEDs are based on observations for large apertures 
taken from the literature. The key rest frame mid-IR range uses the ISOCAM-CVF low 
resolution spectra of Forster Schreiber et al. (2003) for M82, O. Laurent (priv. comm) for 
Arp 220, and ISOPHOT-S (Spoon et al. 2001) and groundbased spectra (Roche et al. 1986) 
for NGC 4418. The colors of the nuclear region of the Type 2 AGN NGC 1068 shifted to 
z~1.67 are indicated by a cross. 
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Fig. 4. — Ratio of 1.2mm and 24/im flux density as a function of redshift for local SED 
templates. The SEDs of M82 and Arp220 represent typical and extremely obscured star- 
forming infrared galaxies. NGC 4418 is an extremely obscured compact object possibly 
hosting an AGN. PDS456 is a luminous quasar whose infrared SED is dominated by the 
AGN mid-IR emission, with a weak far-IR component. The AGN dominated nuclear region 
of NGC 1068 is only shown for a single redshift where the rest frame submm emission of 
the nucleus can be spatially isolated. The horizontal continuous line shows the mean ratio 
for our sample. Horizontal dotted lines indicate the range spanned by the targets that are 
individually detected by MAMBO. The horizontal dashed line indicates the median 1.2mm 
to 24/im flux ratio for submillimeter galaxies (SMGs) - see text for references. 



